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Inspired by a recently proposed model of millicharged atomic dark matter (MADM), we analyze 
several classes of light dark matter models with respect to CoGeNT modulated and unmodulated 
data, and constraints from CDMS, XENON10 and XENON100. After removing the surface con- 
taminated events from the original CoGeNT data set, we find an acceptable fit to all these data 
(but with the modulating part of the signal making a statistically small contribution), using some- 
what relaxed assumptions about the response of the null experiments at low recoil energies, and 
postulating an unknown modulating background in the CoGeNT data at recoil energies above 1.5 
keVee. We compare the fits of MADM — an example of inelastic magnetic dark matter — to those 
of standard elastically and inelastically scattering light WIMPs (eDM and iDM). The iDM model 
gives the best fit, with MADM close behind. The dark matter interpretation of the DAMA annual 
modulation cannot be made compatible with these results however. We find that the inclusion of 
a tidal debris component in the dark matter phase space distribution improves the fits or helps to 
relieve tension with XENON constraints. 



1. INTRODUCTION 

It is an exciting time, albeit puzzling, for direct de- 
tection of dark matter. The three experiments DAMA 
PQ, CoGeNT E] and CRESST @] see hints of light 
DM of mass ~ 10 GeV and cross section ~ 10~ 40 cm 2 
for scattering on nucleons. Although much effort has 
been made to determine whether DAMA and CoGeNT 
best-fit regions (and possibly those of CRESST) can be 
compatible [5]- [21], it is difficult to find any dark mat- 
ter model which achieves this quantitatively, and which 
moreover is consistent with the null results of XENON10 
[23, XENON100 [23, CDMS [2H23, and most recently 
EDELWEISS [2I|. 1 It seems unlikely that all of the 
hints of positive detections can be correct since even non- 
standard DM models with many adjustable parameters 
give a poor fit to the global data [2"5U30| . Arguably, the 
best hope for reconciling the positive and negative signals 
centers around the CoGeNT results, which can be com- 
patible with the lightest and most weakly coupled DM 
models, hence those having the greatest chance of being 
marginally compatible with the null results. This is our 
motivation for focusing the present study on CoGeNT. 

Another of the difficulties with CoGeNT is that it 
presents an annual modulation amplitude that seems to 
be too large compared to the unmodulated rate unless 
a nonstandard velocity distribution is invoked for dark 
matter in the galactic halo [3"lTf3"3"j . This problem is only 
exacerbated by recent results indicating that the ear- 
lier published events included contamination from surface 
events [34 . If CoGeNT is indeed seeing dark matter, it 
is tempting to speculate that there is some contamina- 
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1 Ref. |27l appeared as we were finishing this work, and is therefore 
not included in our analysis. Since it obtains weaker limits than 
XENON, we do not expect this to affect our conclusions. 



tion of the modulating part of the signal that should be 
subtracted, to ameliorate this tension. We explore the 
effect of including such a background, finding that the 
CoGeNT-allowed regions of parameter space are hardly 
affected by the modulating part of the data, regardless of 
any such background, due to their lower statistical sig- 
nificance relative to the time-averaged part of the signal. 

Although we also consider standard elastic and inelas- 
tic dark matter in the present study, one of our motiva- 
tions for carrying it out was the proposal of a new class 
of millicharged atomic dark matter (MADM) models in 
[35] (a simpler variant of models previously discussed in 
refs. [SHI [37]), with emphasis on the special case where 
the dark "electron" and "proton" constituents have equal 
masses, m e = m p . In the ref. [35) . rough estimates were 
made as to the parameter values required to explain the 
CoGeNT observations. Here we better quantify those 
by direct comparison to the data. We will show that 
MADM gives an acceptable fit to the data for the dark 
H atom mass toh ~ 10 GeV, the hyperfine mass splitting 
SrriH ~ 25 keV, and the fractional charge of the e and p 
constituents e ~ 0.01, which controls the scattering cross 
section. The fit to the CoGeNT data by MADM is better 
than the one of standard elastic dark matter (eDM) , and 
only slightly worse than that of standard inelastic dark 
matter (iDM), which requires similar values for the mass 
and mass splitting. 

It is interesting to note that in the special m e = m p 
case of interest, the effective charge of dark atoms is zero 
even at short wavelengths, due to the cancellation be- 
tween the charge clouds of the two constituents, and the 
model becomes an example of magnetic inelastic dark 
matter [55]- [15] (this identification has not been pre- 
viously realized). Magnetically interacting dark mat- 
ter has recently been promoted as a model that can fit 
"everything," all the null and positive direct detection 
constraints @5], including DAMA and CRESST. In the 
present work, we will not reach such optimistic conclu- 
sions, even though we find some room for CoGeNT to 
be compatible with the Xenon and CDMS constraints, if 
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we exercise a moderate degree of skepticism concerning 
the claimed sensitivity of the Xenon experiments to low 
recoil energy events. 2 

We also consider the effect of astrophysical uncertain- 
ties, including nonMaxwellian velocity distributions [50] 
and contributions from debris flow [5TJ|52]. We find quan- 
titative changes to the allowed regions as a result of these 
variations, and in the case of debris flow, some improve- 
ment in the quality of the fits and in the relative size of 
CoGeNT-allowed versus excluded regions. In the remain- 
der of the paper, we will review the theoretical framework 
(section [2]) and our methods of analysis of the various ex- 
periments (section [3]). Results are presented in section [4] 
and conclusions in section [U 



2. THE MODELS AND THEIR PREDICTIONS 

We briefly review the MADM model and highlight its 
differences relative to standard (inelastic) dark matter 
for direct detection. The MADM model contains a dark 
"electron" e and "proton" p bound into dark atoms H 
by a new U(l) gauge force of strength a', whose mass- 
less vector boson kinetically mixes with the photon, re- 
sulting in a small electric charge ±ee for the e, p con- 
stituents. MADM therefore interacts electromagnetically 
with protons in ordinary nuclei. In the special case where 
m e = m 
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the electric coupling vanishes in the leading 
Born approximation and the dominant transition is via 
the atomic hyperfine interaction. It is inelastic due to 
the hyperfine splitting 8mn = |a:' 4 m 2 m 2 /(r7i e + m p ) 3 
which can naturally be of order 10 keV. 

In general, the differential rate of DM scattering on 
nuclei with nuclear recoil energy Er is given by 
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where Nj- is the number of target nuclei, p & = 0.43 
GeV/cm 3 [53] is the DM density in the solar neighbor- 
hood, rriDM is the DM mass, f(v) is its velocity distri- 
bution (see appendix [A] for details), and v m [ n is the mini- 
mum DM velocity needed to produce the given recoil en- 
ergy. In terms of the momentum transfer q — ^JItyi^Er 
for nuclear mass tun and DM-nucleus reduced mass /i, it 
is Wmin = ?/ (2fi) + 8m/q for inelastic DM with mass split- 
ting 5m. The differential cross section for DM of velocity 
v scattering on nuclei can be expressed as 



der 



m N a p 
2^/4 



z 



In 

Jv 



(A-Z) 



F 2 N (q)F 2 DM (q,v) (2) 



where a p is the cross section for DM scattering on protons 
at threshold, jj, n the DM-nucleon reduced mass, F/v the 



2 We are not willing to consider 7a-8a allowed regions of DAMA 
for the purpose of finding overlap, as done in |49l 



Helm nuclear form factor [53], and Fum the DM form 
factor. For generic DM, f n /fp = Fr,m = 1, whereas for 
MADM, /„ = and F DM ~ (v 2 - v 2 min )/q 2 . 

In the MADM model, the differential detection rate 
can be expressed as 
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where the velocity integral is given by 
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Because of the Fum form factor, a p cannot be defined 
in the usual way as the cross section at threshold, but if 
desired, one can define an effective a p by inserting char- 
acteristic values vq, qo of the DM velocity and momentum 
transfer, 
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for example go = 25 MeV (appropriate for scattering on 
germanium at Er = 5 keV), vq = 220 km/s. The form 
factor has compensating such terms, 
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so that the rate is independent of ^o, qo- 

In the remainder of the paper, we will be interested 
not only in MADM but also in generic DM models where 
Fdm = 1, which can be elastic (f m in = 9/(2/^)) or in- 
elastic (w min = q/(2fi) + Sm/q). We note that MADM 
differs from the generic inelastic DM model (iDM) not 
only because of the form factor but also because of 
its mild isospin violation, f n — 0. In this work we do 
not investigate models with arbitrary isospin violation, 
such as the case of /„/ f p = —0.7 that significantly weak- 
ens the sensitivity of xenon-based experiments [HI [28]- 
I5U1 155H5T] . However we will consider both the isospin 
conserving eDM and iDM models, and the special case 
in which /„ = 0. The latter occurs naturally when DM 
interactions are mediated by a U(l) gauge boson that 
kinetically mixes with the photon. 



3. METHODOLOGY 

In this section we describe the computation of x 2 f° r 
the CoGeNT modulated and unmodulated data, DAMA 
modulations, and the exclusion contours coming from 
XENON100, XENONI0 and CDMS-Ge (CDMS-Si gives 
weaker limits for the models of interest). Some of this 
is standard, following the same procedures as the rele- 
vant experimental collaborations themselves. In the case 
of CoGeNT, we suggest an additional background to be 
subtracted (over and above the surface contamination 
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events discussed in 34J), namely an unspecified modulat- 
ing background that we take to be independent of energy, 
to fit away the modulation of the CoGeNT high-energy 
bin f .5-3.1 keVee. However in the end we will find that 
the experimental errors on the modulating part of the 
signal are too large for any of these details to matter 
for the combined fit to the modulating and unmodulated 
data. 



3.1. CoGeNT unmodulated signal 

We start by fitting to the CoGeNT unmodulated signal 
by itself, since this has higher statistical significance than 
the annual modulation signal. We analyzed the public 
release data [SB] provided by the CoGeNT collaboration 
by reweighting each event using the detector efficiency 
fefi(E) [3]. To remove the known background due to cos- 
mogenic L-shell electron capture events in the region of 
interest, we used the fitting parameters from the K-shell 
peaks provided by CoGeNT [58 and the ratio between L- 
shell and K-shell decays [59] (see e.g. also [30l|3T] for the 
parameters of the cosmogenic backgrounds). To relate 
the observed ionization energy to the inferred nuclear re- 
coil energy, we take the same quenching factor as used by 
CoGeNT [58], E/keVee = Q(E R /keVnr)P , with Q = 0.2 
and p = 1.12. 

Subsequent to their published results, CoGeNT re- 
ported a new estimate of the surface event contamina- 
tion near the energy threshold [33] . We apply the central 
curve of the estimated surface event correction factor on 
p. 14 of ref. |34j to the energy spectrum after removing 
the cosmogenic background. A constant background of 
2.685 cpd/keVee, taken as the average of the event rate 
between 2 keVee and 3 keVee, is further subtracted in 
order to remove the events above E = 1.5 keVee, whose 
energy dependence is approximately constant. This is 
motivated by the fact that light DM models that account 
for the large excess at lower energies predict essentially 
no events in this higher energy range. For the errors on 
each bin, we consider the statistical errors based on the 
original events in the bin before efficiency correction and 
background subtraction, and also the uncertainties in the 
parameters of the cosmogenic background. 

In fig. [T] we display the data with the cosmogenic back- 
ground, surface events contamination, and the constant 
background subtracted, along with the prediction of the 
MADM model with parameter values ibh = 9.9 GeV, 
S u = 24.7 keV, e = 0.029. Below we will show that this 
model (which we refer to as our benchmark) corresponds 
to the best fit parameters and is marginally allowed by 
XENON constraints, as we will determine them. 

To accurately predict the theoretical DM contribution, 
one should take into account that the 442 day period of 
CoGeNT data-taking is not an integer multiple of a year, 
as well as account for inspection and power outage days of 
the CoGeNT operation, and average the modulated rate 
over the actual period rather than assuming the mean 
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Figure 1: CoGeNT energy spectrum (subtracting cosmo- 
genic background, surface contamination, and constant back- 
ground), and the prediction from the MADM model, (?tih, 
<5h, e)=(9.9 GeV, 24.7 keV, 0.029). 
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Figure 2: CoGeNT data and MADM (the same model as in 
Fig.|l| prediction for modulation in the low-energy 0.5-1.5 
keVee bin. Here we also consider the possibility of including 
a modulation background (dashed line) on top of the MADM 
signal in order to explain the CoGeNT modulation events in 
the high energy bins. The modulation background used here 
is Dsin((27r/year)(i - to)) where D = -0.487 cpd/kg/keVee 
and to = 199 d. As shown in the figure, the CoGeNT modu- 
lation amplitudes are larger than MADM predictions. 



value of the earth's velocity relative to the DM halo rest 
frame. We find that the latter makes a —1.7% correction 
in the predicted rate, too small to be of interest given the 
current experimental uncertainties. 



3.2. CoGeNT annual modulation 

We analyzed the modulated data obtained from the 
CoGeNT collaboration, following ref. [31] by grouping 
the events into a low-energy bin with E < 1.5 keVee and 
a high-energy bin between 1.5 and 3.1 keVee. To deter- 
mine the modulation spectrum, we model the cosmogenic 
background by taking into account the outage days. We 
also remove the constant background and the estimated 
surface contaminated events |34j assuming that they are 
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Figure 3: (a) Left: comparison of 90% c.l. limits on elastic DM cross section derived by XENON100 [23] using profile likelihood 
method (solid, red curve) and by us using maximum gap method with central value (dashed, blue curve) and 2a lower limit 
(dot-dashed, black curve) of the XENON100 relative scintillation efficiency C e g. (b) Right: XENONfOO acceptance curves 
A(Er) for the central values (solid lines) and 2a or la lower limits (dashed lines) of £ ff, taken from ref. [23] (left curves) and 
|60| (right curves). The DM mass is taken to be 10 GeV for the choice of cut acceptance. 



unmodulated in time. This is supported by the fact that 
the surface events originally rejected by the CoGeNT 
analysis of ref. [3] do not show any sign of modulation. 
Based on the fit to the unmodulated data, where no dark 
matter events were found to occur above 1.5 keVee, we 
consider the modulating signal in the high-energy bin to 
be due to an unspecified background. To model it using 
the smallest number of parameters, we take this extra 
background to be independent of energy, having a rate 
of the form 

D sw(u)(t - t )) (7) 

Assuming that ui — 2ir y _1 so that only D and the phase 
to are free parameters, we find that the best fit values 
are D = -0.487 cpd/kg/keVee and i = 199 d. It will 
be seen that these values do not change very much when 
we allow them vary in the minimization of the total \ 2 
for all the CoGeNT data. The energy dependence of this 
hypothetical background is unknown, so we considered 
the two cases where it either reduces only the signal in 
the high-energy bin, or that in both. We did not see 
significant differences between the two possibilities in the 
final results. 

In any case, one needs to separate the data into an 
unmodulated part (which has already been considered 
above) and the modulating component whose contribu- 
tion to x 2 is computed separately from that of the un- 
modulated signal. We compared several different ways 
of doing this: fitting the time dependence in a given 
bin to a constant plus sinusoidal term; splitting the data 
into its time-averaged value plus fluctuations; or defining 
the time-dependent part by subtracting the time average 
only over consecutive 365 day periods rather than 442 
days (which should be more correct). However, the data 



are not yet at the level of precision where these distinc- 
tions are important. We find that the time-averaged part 
of the rate in the low-energy bin is 2.47 ± 0.1 cpd/kcVce 
(corresponding to the range Sx 2 = 1)- 

The raw CoGeNT modulation data from the low- 
energy bin are plotted in fig. [2] along with the MADM 
prediction plus modulating background |7| for the same 
model parameter values as used in fig. II] Although the 
amplitude of the theoretical prediction looks small com- 
pared to the central values of the data, the experimen- 
tal error bars are large, and so the contribution to the 
total x 2 horn this discrepancy is relatively small. For 
the benchmark model, (mn, 8n, e)=(9.9 GeV, 24.7 keV, 
0.029), removing the expected DM contribution from the 
modulation signal in the low-energy bin (0.5-1.5) keVee 
(i.e., keeping only the modulating background ^ as sig- 
nal in this energy range), corresponds to a reduction in 
the total CoGeNT A X 2 of only - 0.5. 

To display the regions of parameter space preferred 
by the data (see section ffl), we use the A^ 2 = 
3.53, 8.03, 14.2 contours, corresponding to confidence 
intervals of 68.3%, 95.5%, and 99.7%, appropriate for fit- 
ting three model parameters (ran, Su, e). 



3.3. XENON100 

The XENON100 collaboration [23] observes 3 events 
at nuclear recoil energies of 12.1, 30.2 and 34.6 keV re- 
spectively, all of which are too high to be directly ac- 
counted for by our best-fit models. We compute the num- 
ber of events above the XENON100 low-energy thresh- 
old 8.4 keVnr using Poisson fluctuations [ST]: N — 
J dE R dR/dE R A(E R ) times the exposure 100.9 kg x 48 
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Figure 4: Solid curves: limits from XENON100 assuming the the central and 2er lower limits of the XENONfOO £ ff (leftmost, 
red lines) and central and lcr lower limit of the ZEPLIN-III £ e fl (rightmost, blue lines). Dashed curves (black): limits from 
XENON10 assuming the three possible electron yields versus energy shown in fig. [5] The shaded areas are the CoGeNT best-fit 
regions for eDM (left graph) iDM (central graph) and MADM (right graph) models, discussed in section El 



d. The acceptance function A{Er) depends upon data 
quality cuts and S2/S1 energy resolution discrimination 
cuts which we take from [23]. For the data quality cut 
acceptance, a linear interpolation based on XENONlOO's 
analysis of the DM mass dependence is used. We elabo- 
rate upon this below. 

Ref. [23] uses the profile likelihood method [52] to de- 
rive its limits on dark matter. It would be difficult to 
reproduce these results without having full access to the 
data (including backgrounds). As a simpler alternative, 
we adopt the maximum gap method of ref. [63] . In this 
method, the energy range of the search window, 8.4-44.6 
keVnr, is partitioned into 4 intervals ("gaps") separated 
by the three observed events: 8.4-12.1, 12.1-30.2, 30.2- 
34.6, 34.6-44.6. One then looks for the interval in which 
the observation of no events is the least probable for a 
given DM model, and uses this to set a limit on the cross 
section. Not surprisingly, for the low-mass DM models 
we are interested in, this turns out to be the lowest in- 
terval, 8.4-12.1 keV. 

We determine the expected number of events in each 
gap including not only the dark matter signal but also the 
expected background, which is estimated to be 1.8±0.6 
events 23 in the full window. We assume this back- 
ground is uniformly distributed in energy. Since nearly 
all the light DM events in the full window are within the 
first gap, the DM contribution to the number of events 
x expected within this gap is nearly the same as the DM 
contribution to the total number /i in the full search 
window, and the function Co(x,/j,) that gives the Pois- 
sonian probability that more than events should have 
been observed in this window simplifies to Co(x,fi) = 
1 — e~ x (/z — x + 1). The 90% c.l. limit is determined by 
setting Co = 0.9. Both x and /i get contributions from 
DM and from the background, but the DM contribution 
to each is nearly equal for light DM. Therefore \i — x is de- 
termined by the number of expected background events 
outside the first gap, ijl — x = 1.62. The number of DM 



events allowed at the 100 Co% confidence level is then 
given by x DM = ln(2.62/(l - Co)) - 0.18, where 0.18 is 
the number of background events expected in the lowest 
gap. For Co = 0.9 this gives xdm =3.1. 

In fig. J3ja) we show that the maximum gap method 
gives reasonable agreement with the more sophisticated 
profile likelihood method at masses above 8 GeV. At low 
masses, the maximum gap method gives a stronger limit 
than the XENON100 result [23 , unless wc weaken the 
assumed sensitivity by adopting the 2er lower limit of 
XENONlOO's relative scintillation efficiency £ e ff (more 
about which below) . This is not surprising since the pro- 
file likelihood method essentially marginalizes over £ e ff, 
which gives a better fit to light DM models (weaker limit) 
for smaller £ e ff ■ In the higher mass range, the maximum 
gap method gives a limit that is close to the one derived 
from XENONlOO's central £ e ff- We will need to weaken 
this limit as described presently in order to accommodate 
the CoGeNT-preferred region. But it should be kept in 
mind that the additional weakening we invoke is not so 
great a departure from XENONlOO's central £ e ff value 
as one might at first think, for dark matter masses below 
9 GeV. This is because of our finding that the maximum 
gap method yields stronger limits in this mass range, for 
a given choice of C c g, than does XENONlOO's profile 
likelihood approach. 

To compute the expected number of DM events in a 
given energy interval, one needs to weight the raw event 
rate dRj 'dE(Efi) by the acceptance function A(Er), de- 
tails of whose computation are given in appendix [B] The 
acceptance at low recoil energies depends strongly upon 
the relative scintillation efficiency £ e ff at low En [64] [65] . 
The C cS {E R ) used by XENON100 in their own analysis 
is measured down to E = 3 keV [6B] and logarithmically 
extrapolated to zero for lower energies. Previous mea- 
surements of £ c h(Er) had only gone as low as 4 keV. 
The validity of the new measurement at the lowest en- 
ergy has been questioned in ref. |67) . which claims that 
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it should be regarded only as an upper limit on the true 
efficiency of the XENON100 detector. 

It will be seen that XENON100 excludes most of the 
parameter space preferred for fitting the CoGeNT data 
unless some modest reduction of the assumed sensitivity 
at low recoil energies is made. In the spirit of our "opti- 
mistic" analysis, and taking into account the reservations 
mentioned above, we adopt the recent measurement of 
£efi by the ZEPLIN-III collaboration [50], which is also 
used in ref. [57], and supersedes the earlier ZEPLIN-III 
measurement |68] that was criticized in ref. [641 . This 
determination of £ e ff goes down to Er — 2.2 keV. We 
follow the practice of ref. [H] in linearly extrapolating it 
from this point to zero at En = 1 keV to give some rea- 
sonable estimate of the response at lower energies. Using 
an £ e ff close to this one, we will show that substantial 
overlap between the XENON100 and CoGeNT allowed 
regions can be achieved. In section [4] we will display the 
XENON100 constraints using the 0.5a lower limit (de- 
rived from 60J) of the measured £ e g(E). To illustrate 
the impact these different assumptions make on the sensi- 
tivity of the XENON100 experiment to low-recoil events, 
we plot in fig. J3jb) the full acceptance function at low 
recoil energies, based upon the £ e ff used by XENON100 
[2"3] (along with its 2a lower limit) and that of ZEPLIN- 
III (showing also its la lower limit). As we will discuss 
below, our benchmark MADM model (ran, <5h, e)=(9.9 
GeV, 24.7 keV, 0.029) predicts only 2.8 events using the 
ZEPLIN-III 0.5a lower limit C e g, which is less than the 
XENON100 90% c.l. upper limit of 3.1 events. 

To show the effect of the different possibilities for £ e ff 
on the allowed DM parameters, we display in figure [4] 
the resulting XENON100 limits on the eDM, iDM and 
MADM models for the central and 2a lower limits of 
the XENON100 C c g (solid red curves), as well as for 
the central and la lower limit of the ZEPLIN-III C c s 
(solid blue curves). The ZEPLIN-III la error bars are 
comparable in size to the 2a range of XENON100 since 
the latter is based upon an average over several different 
experiments. For reference, the CoGeNT-allowed regions 
for eDM, iDM and MADM models are also shown (these 
will be discussed in section [4| . Although the CoGeNT 
regions are ruled out by the XENON100 central and 2a 
C e gs, and and marginally by the ZEPLIN-III central C c g, 
one sees that there is large overlap with the ZEPLIN-III 
la C cS , and even with the ZEPLIN-III 0.5a C cS . To be 
conservative in our optimism, we will adopt the ZEPLIN- 
III 0.5a £ e ff for the remainder of our analysis, but it 
should be kept in mind that one could reasonably argue 
for a case that makes even greater room for a viable DM 
explanation for the CoGeNT events. 



3.4. XENON10 

Unlike XENON100, the XENON10 experiment relin- 
quishes its demand for scintillation photoelectrons (SI) 
to veto electron recoils, in order to gain a lower energy 
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Figure 5: Different possibilities for the relation between nu- 
clear recoil energy and number of S2 electrons observed by 
XENON10. Top curve (solid) is from ref. [B7J EH], bottom 
curve (dotted) is that used by [22] . The intermediate (dashed) 
curve is our hypothetical ansatz that allows for compatibility 
with the CoGeNT events. 



threshold and thus greater sensitivity to light dark mat- 
ter. Thus one does not require a significant acceptance 
function, and need only consider the raw rate of nuclear 
recoils in energy intervals bounded by events within the 
range 52 = (5,43). This corresponds to recoil energies 
En in the range (-Emm, E max ) where E mm is associated 
with the minimum number of 5 electrons imposed by the 
experiment as a cut on the data, and E maK = 10 keV, 
whose exact value is not important for constraining light 
DM. XENON10 observes 23 events in this range of the 
S2 signal. 

To set a limit, energies are assigned to each event ac- 
cording to the relation S2 = EnQ y (Ejf), where Q y (En) 
is the electron yield function that must be estimated from 
a theoretical model. Then expected numbers of events 
Xi in every energy interval bounded by two events are 
compared to measured events Ni in these intervals (ex- 
cluding the contributions from the endpoints), allowing 
for the computation of a probability pi that more events 
should have been seen than measured, based upon Pois- 
son statistics. Finding the maximum of pi over all in- 
tervals and comparing to the known function pMax(C, /i) 
leads to excluded regions in the parameter space at a 
specified confidence level C, where fj, is the expected num- 
ber of events in the interval (-E mm , -Emax)- This is the 
Pm&x method of ref. [63], which was used by XENON10 
to set limits on light dark matter in [52] , and in the anal- 
ysis of ref. [30]. We too use this method to compute the 
90% c.l. limits from XENON10. 

The crucial question is what the actual form of the elec- 
tron response function Q y should be, and the resulting 
threshold energy E m j n . XENON10 relies upon an extrap- 
olation of Lindhard's theoretical prediction below 4 keV, 
which leads to the estimate E m i n = 1.4 keV. This has 
been criticized by refs. [67j[69], which advocate a differ- 
ent prediction that leads to lower estimated sensitivity 
of the experiment to light DM events, and the higher 
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threshold -Emm = 3.6 keV. We illustrate the differences 
between the two estimates in fig. [5] The XENON10 ver- 
sion of Q y leads to a strong conflict with the regions 
of parameter space preferred for explaining the CoGeNT 
signal, while the relaxed assumption removes this tension 
entirely. We will consider a hypothetical intermediate 
ansatz, also shown in fig. [5] that falls between the two 
estimates. In fig. [4] we show the XENON10 limits that 
result from the three choices of Q y corresponding to fig. 
[5j these are the dashed (black) curves. Comparing to the 
CoGeNT-allowed regions, one sees that it is sufficient to 
choose the hypothetical intermediate Q y ansatz for main- 
taining compatibility between CoGeNT and XENON10. 
We adopt this choice in what follows. 



3.5. CDMS 

For CDMS-Ge, we follow the procedure of ref. [3U] in 
assigning a x 2 of the form 

36 

where N is the number of DM events in the relevant en- 
ergy range, from the low-energy threshold of 2 keV up to 
20 keV, using the efficiency for the best-performing detec- 
tor T1Z5 as provided by the experiment. The ability of 
CDMS to distinguish 2 keV recoils from background has 
been disputed in ref. [70] , which argues that 5 keV could 
be considered as a realistic threshold considering the un- 
certainties. However we are not forced to consider such 
a weakening of the ensuing CDMS constraints since even 
with 2 keV threshold, it will be seen that they do not con- 
flict with the CoGeNT preferred regions. The previously 
observed tension between CoGeNT and CDMS-Ge has 
been alleviated by the subtraction of the CoGeNT sur- 
face contamination events, which has reduced the former 
"signal" events by ~ 70% near the low energy thresh- 
old. In addition to including ^ in our total x 2 (al- 
though it makes no contribution in the relevant regions) , 
we plot the 90% upper limit from CDMS-Ge correspond- 
ing to 45.9 events. We have also analyzed constraints 
from CDMS-Si data and found them to be weaker, and 
therefore do not consider them henceforth. 

Recently CDMS has also placed limits on a modulating 
DM signal in the energy range 5-12 keVnr [26]. However, 
considering the CoGeNT quenching factor, this essen- 
tially coincides with the high-energy bin in which our 
preferred models predict no events. Therefore we are not 
constrained by these data, and do not consider them in 
our analysis. 

We note that an independent analysis of CDMS data 
by ref. [7T] finds an excess of low-energy events that is 
consistent with the CoGeNT signal. This finding is also 
consistent with our own observation that the CDMS con- 
straint is the weakest of the null results, and that we 
did not have to make any special assumptions to achieve 
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9.36 
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9.47 


total 


70.9 
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Table I: The best-fit models for CoGeNT in standard eDM, 
iDM, and in MADM models. All the models in the table 
have isospin violation /„ = 0. The case of /„ = f p for eDM 
or iDM is the same except for decreasing a xv by a factor of 
(A/Z) 2 = 5.15 = 10 on . 



compatibility with the CoGeNT allowed regions that we 
determine below. 



3.6. DAMA 

We have computed the effect of DAMA's annual modu- 
lation data on our fits, following the same procedure as in 
ref. [30]. Unlike the CoGeNT allowed region, which can 
be made compatible with even the most stringent con- 
flicting data by making the judicious choices discussed 
above, the DAMA allowed regions are firmly ruled out. 
The minimum value of x 2 = 68 without DAMA increases 
to 152 when DAMA data are included, even though only 
8 data points (the number of energy bins) are included. 
For this reason we do not try to reconcile both CoGeNT 
and DAMA results with the Xenon constraints. We con- 
clude that the DAMA signal is due to some unknown 
background in the context of our models, as other recent 
studies also seem to indicate [72]. 

4. RESULTS 

4.1. Standard DM velocity distribution 

We determined the regions of parameter space that 
give the best fits to the three models under considera- 
tion: standard elastic and inelastic dark matter (denoted 
as eDM and iDM, respectively), and the millicharged 
atomic (MADM) model. For eDM, there are just two 
model parameters, the DM mass m x and cross section 
on nucleons a n , while iDM also includes the mass split- 
ting 8 X . For MADM, we use the fractional charge e to 
parametrize the interaction strength instead of <r n , and 
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Figure 6: CoGeNT best-fit regions (shaded) in the m x -a xp 
plane, and "optimistic" 90% c.l. upper limits from XENONfO, 
XENONfOO and CDMS-Ge for eDM (top row) and iDM 
(bottom row) models, with isospin violation /„ = (left- 
most graphs) or isospin conservation (rightmost graphs). The 
shaded regions correspond to 68.3%, 95.5% and 99.7% confi- 
dence intervals. 



we distinguish the mass and mass splitting from those of 
standard DM using mn and dmn (where H stands for 
the dark "hydrogen" atom). For all models, we adopt 
the relaxed versions of the XENON100 and XENON10 
constraints corresponding to the ZEPLIN-III 0.5-<r lower 
limit on £ c ff (fig. |4| and the hypothetical intermediate 
Q y (fig. [5]), respectively. 

In table |T] we display the parameters corresponding to 
the models that minimize the \ 2 f° r the CoGeNT data 
(as well as CDMS although the latter gives a vanishing 
contribution to \ 2 at these best fit points), as well as the 
corresponding minimum values of \ 2 - Inelastic DM with 
m x = 9 GeV and S x ~ 24 keV gives a significantly better 
fit than does elastic, and the millicharged atomic dark 
matter model, also inelastic, gives nearly as good a fit 
as iDM, with a slightly larger mass mn — 10 GeV and 
mass splitting <5h ~ 25 keV. Due to the relatively large 
inelasticity, the DM-nucleon cross section in the iDM case 
is about 3 orders of magnitude higher than that of eDM. 

For the eDM and iDM models, fig. [6] displays the 
CoGeNT 68.3%, 95.5% and 99.7% confidence intervals, 
along with the constraints from XENON10, XENON100 
and CDMS, in the plane of m x and the cross section 
for DM scattering on protons. The XENON100 limits 
are the most constraining, cutting through the center of 
the CoGeNT allowed regions, but still leaving room for 
compatibility within the 68.3% c.l. intervals. The best-fit 
points are excluded by XENON100 for eDM but allowed 
for iDM. Having untuned isospin violation of the form 



/„ = (DM coupling only to protons) slightly weakens 
the Xenon constraints relative to the CoGeNT regions. 

Next we turn to the MADM model, which we study 
in somewhat greater detail. The best-fit parameter set 
for this model is marginally consistent with XENON100, 
and we adopt it as our benchmark model. In table [TT| we 
present the contributions for this model to \ 2 from the 
nonmodulated and modulated data and from comparison 
to the DAMA annual modulation, and numbers of events 
expected for XENON and CDMS experiments. In figure 
[7j the best-fit regions and 90% c.l. constraints for MADM 
are displayed in the three planes mu-e, Su-e, and Jh-wh, 
where in each case the parameter that is not varied is 
fixed to its benchmark value. Approximately half of the 
CoGeNT 68% confidence intervals are left unexcluded by 
the XENON100 limit, which is the most stringent of the 
constraints. 
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Table II: Contributions to the total \ 2 or to the number of 
events observed by CDMS-Ge or XENON experiments, for 
the benchmark MADM model with (mn, 6n, e)=(9.9 GeV, 
24.7 keV, 0.029). For the CoGeNT modulation spectrum, we 
include the modulating background {7j|. For the XENON100 
and XENON10 limits, we consider various choices for £ e ff and 
Q y , discussed in section [3] 



4.2. Nonstandard DM halos 

We have investigated the effects of two possible modi- 
fications to the standard Maxwellian phase space density 
of the dark matter. The first is a distortion of the s hape , 



of the type f(v 2 ;v 2 csc ) -> f k {v 2 /k; v 2 c Jk) (see eq. 
that has been suggested in ref. [50] based upon Edding- 
ton's relation between f(v) and and the density profile 
p(r). Comparisons to V-body simulations and to the 
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Figure 7: CoGeNT best-fit regions (shaded) for the MADM model near the benchmark point (m H , 5h, e)=(9.9 GeV, 24.7 keV, 
0.029). The shaded regions correspond to 68.3%, 95.5% and 99.7% confidence intervals. In each plane, the parameter that is 
not varying is fixed to its best-fit value. 




Figure 8: Dependence of best fit values of mH, Sn, e and 
minimum value of x 2 on the halo velocity profile index k, in 
the MADM model. 



Milky Way halo give k — [1.5,3.5] and k ~ 2 respec- 
tively. The effects of varying k on fitting to CoGeNT 
were previously studied in ref. |31j . We find that increas- 
ing k has a relatively small effect on the best-fit regions 
and on the constraint curves, shifting them to somewhat 
higher values of DM mass and cross section, and smaller 
values of the mass splitting. We illustrate this for the 
MADM model in figure [HJ where the best-fit values of 
the parameters and the minimum value of \ 2 are plotted 
as functions of k. 

We next consider the debris flow scenario of ref. [52] • 
Tidal debris is a form of unvirialized dark matter due 
to tidal stripping of ubiquitous subhalos in the galaxy. 
Unlike tidal streams, which are rare features whose ef- 
fects have been previously considered for interpreting Co- 
GeNT in ref. [72], debris flow is argued in [52 j to be 
homogeneously distributed in the galaxy and therefore 



guaranteed to be present in the solar system, with a dis- 
tribution that was inferred in |52] by analysis of A^-body 
simulation results. 

We correct an error [73] in [52] by writing the full dis- 
tribution function as 



f(v) = (l-e )f v (v) + e f d (v) 



(9) 



where /„ is the Maxwellian distribution of the virialized 
component, fd is the contribution from debris, and eo = 
0.22 is determined by comparison to A-body simulations. 
In |52j . the coefficients of f v 4 were mistaken for functions 
of velocity in the computation of the DM detection rate. 
Ref. [52] further models fd{v) (in the earth frame) by an 
isotropic function with support in the window Vf — v e < 
v <v f + v e , 



h{v) = 



2vf v, 



■Q(vf + v e - v)Q(v - v f + v e ) (10) 



where v e (t) is the modulating earth velocity and Vf is 
a flow velocity characterizing the debris. (f v ,d are both 
normalized such that J dvf v .d = 1.) In [52], the value 
Vf = 340 km/s was suggested, but we find that a higher 
value Vf = 475 km/s is required for consistency with 
their determination of the ratio e(v m - ln ) of debris parti- 
cles to the total number of DM particles having v > v m i n . 
In fig. |9|we show that e(v m i n ) as directly computed from 
) is in reasonable agreement with the result from 
on analysis of A^-body simulation data, for our 
choice oi Vf, which is definitely not the case for Vf = 340 
km/s. This provides a consistency check on the approxi- 
mation (10) with Vf — 475 km/s. The constant eo in ((9l 



igik 

eqs. ( |9jl0 
basec 



is equal to e(0). 

The debris flow contribution has interesting effects on 
the quality of fits and the best-fit regions for dark matter 
direct detection. We summarize the best-fit values in 
table |III| The graphical results are presented in fig. [lOj In 
the case of isospin-conserving elastic DM, the goodness of 
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Figure 9: The ratio of debris particles to total number of 
dark matter particles having v > u m i n in the earth frame, as 
determined in ref. [52] (solid and dot-dashed curves), and as 
calculated here, using the approximation ( |10[ | in eq. (j9|. 
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Table III: The best-fit models for CoGeNT in standard eDM 
and iDM and MADM models, including the debris flow com- 
ponent in the DM phase space distribution function. The 
eDM and iDM models are isospin conserving, f n — f p . Num- 
bers in parentheses refer to the second-best fits in case of a 
second local minimum of \ 2 . 



fit is nearly unchanged, but the best-fit mass is shifted to 
the lower value to# — 7.4 GeV, while the cross section 
increases 3 by a factor of ~ 5, thus moving below the 
region excluded by XENON100, whose constraints are 
also shifted, but by a lesser amount. 

For inelastic DM, debris flow leads to an improve- 
ment in the best fit with \ 2 = 59.8 (smaller by 3.6 than 
the Maxwellian halo result). The most dramatic change 
amongst the best-fit model parameters relative to the 
Maxwellian halo values is in the lower value of the cross 
section. There is no noticeable relative shift between the 
CoGeNT best-fit region and the Xenon constraints in this 
case. However a new local minimum with x 2 = 63.3 de- 
velops at lower m x , higher S x and larger a xp . It origi- 
nates from the presence of two components in /, only one 
of which dominates the signal near a given local mini- 



3 Recall that one must divide the best-fit u X:P values in table [j] by 
5.15 to convert them to the isospin-conserving case 



mum. The new low-mass region requires higher recoil 
velocity for a given recoil energy (especially at large mass 
splittings) , and the debris flow contribution dominates at 
high v. A hint of this second minimum is visible in the 
lower middle graph of fig. 10 which has S x = 20.7 keV, 
the value corresponding tothe global minimum, rather 
than the S x of the secondary minimum. It is clearly vis- 
ible in the bottom middle graph where S x = 24 keV. An 
interesting feature of the second minimum is that, even 
though its x 2 is higher by 3.5 than at the global min- 
imum, it evades the Xenon constraints somewhat more 
robustly. 

The MADM model with debris flow shares qualitative 
features similar to those of iDM. The best-fit \ 2 is better 
(though only by 0.8) than in the Maxwellian case, the 
DM mass is higher, and the fractional charge e is smaller. 
There is still a second local minimum, but in this case 
its mass splitting and \ 2 are closer to that of the global 
minimum, so that both regions appear more clearly on 
the same contour plot at the global 5n best-fit value. 
Again, the secondary minimum at lower mass and higher 
cross section is less constrained by XENON. 



5. CONCLUSIONS 

In this work we have attempted to find the most credi- 
ble set of assumptions that would allow for a dark matter 
interpretation of the CoGeNT excess, without any fine- 
tuning of isospin- violating couplings, while remaining at 
least marginally consistent with constraints from null 
experimental results, mainly those of XENON100 and 
XENON10. Our focus was on standard elastic and inelas- 
tic DM (possibly having mild isospin violation /„ = 0), 
and on a particular example of magnetic inelastic DM, 
the millicharged atomic DM model. Some elements of 
our analysis that are new relative to previous studies are: 
subtraction of recently identified surface-contaminated 
CoGeNT events; hypothesis of a modulating background, 
especially at high recoil energies; and consideration of an 
unvirialized DM debris flow component (correcting an 
error in the literature). 

The main arguments for believing that the DM inter- 
pretation of CoGeNT can still be viable are the experi- 
mental or theoretical uncertainties in the Xenon detector 
response sensitivities (parametrized by the functions £ e ff 
and Q v respectively, for XENON100 and XENON10) at 
low recoil energies. For XENON10 to exclude the Co- 
GeNT preferred region, a theoretical model for Q y is 
required [22] , which is argued to be no more plausible 
than an alternate model [H71 EH] that leaves room for the 
CoGeNT allowed region. We showed that an ansatz for 
Q y intermediate between these two extremes is also ade- 
quate for sufficiently softening the XENON10 constraint. 
Similarly XENONlOO's ability to exclude the CoGeNT 
region relies upon a particular measurement and extrapo- 
lation to low energies of the relative scintillation efficiency 
£ c ff . We showed that the assumption of an alternate £ e ff, 
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Figure 10: CoGeNT best fit regions and constraints for the eDM, iDM and MADM models (left to right) for Maxwellian DM 
velocity distribution (top row) and including debris flow contribution (second and third rows). The middle row sets the mass 
splittings to the values corresponding to the global minima of the \ 2 f° r iDM and MADM, while the bottom row corresponds 
to the higher local minima. The eDM and iDM results are for the isospin-conserving f n — f p case. 
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well within the error bars of another recent measurement 
|60| . is sufficient to make the CoGeNT best- fit model 
marginally compatible with the XENON100 constraint. 

While standard elastically scattering DM can still give 
a reasonable fit to all the data, inelastic models do bet- 
ter. Our featured MADM model has a different depen- 
dence upon the recoil momentum and DM velocity than 
standard iDM, but fits the data nearly as well as stan- 
dard iDM. Our assumption of a hypothetical modulating 
background to explain CoGeNT modulations above 1.5 
keVee improves the % 2 s somewhat, but our results depend 
quite weakly upon this feature, due to the lower statisti- 
cal significance of the modulations compared to the un- 
modulated signal. Distortions from the Maxwellian DM 



distribution functions of the form (A3 1 have a small ef- 
fect on our results. However an additional contribution 
from debris flows (including a correction original to this 
paper) noticeably lowers the DM mass and relaxes the 
tension between eDM and XENON100, improves the fit 
of iDM, and has the interesting effect of introducing a 
second local minimum in the x 2 of iDM and MADM at 
lower DM mass, whose tension with Xenon is reduced 
relative to the global best-fit model. 

Although our results for standard eDM and iDM mod- 
els may be of primary interest for most readers, our orig- 
inal motivation was to study the validity, as an expla- 
nation for CoGeNT, of MADM, which appeared to be 
a new class of model; only in the course of this work 
was its similarity to other inelastic magnetic DM mod- 
els noticed. But the existence of constituents with frac- 
tional electric charge e ~ 0.01 is a distinguishing feature 
of MADM, and a potentially worrisome one since early 
studies of similar models claimed that e > 10~ 3 is ruled 
out [75] by searches for anomalously heavy isotopes of 
hydrogen. In ref. [35] we pointed out several possible 
loopholes to this constraint, including the fact that the 
screening effect of the earth's magnetosphere on the sub- 
dominant ionized component of the dark matter has not 
been taken into account in its derivation. Subsequent 
to the publication of |35j . we have received confirmation 
from some of the authors on the heavy isotope search 
papers (76] [77] that the sensitivity of these searches to 
isotopes with charges deviating from exact integer values 
could have been significantly reduced [75]. We therefore 
keep an open mind about the viability of MADM even 
with fractional charges as large as ~ 0.01 as we find nec- 
essary for explaining CoGeNT. It could be interesting to 
reinvestigate in more detail what values of e are consis- 
tent with existing searches for charged relics. 

Based on these results, we remain optimistic that Co- 
GeNT may really be detecting dark matter of mass ~ 10 
GeV, and perhaps with inelastic couplings between states 
split by ~ 20 keV. New data expected from CoGeNT in 
the near future should shed further light on this exciting 
possibility. 
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Appendix A: dark matter velocity distributions 

The velocity distribution function for the standard 
Maxwellian halo is often taken to be 
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with 



A^ 1 = AttvI 



rf (a) - ae" a - + -a 



1 1 



(Al) 
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where a = v csc /v . Standard values are Vq — 220 km/s, 
v csc = 500 - 600 km/s, v e = 232 ± 30 where ±15.3 is the 
range of annual variation. (The earth's velocity relative 
to the DM barycenter appears in f(v + v e ) when com- 
puting nuclear recoil rates.) Unless otherwise noted, we 
adopt w csc = 600 in computing detection rates. 

In addition to this standard form, we also consider a 
generalized expression with the extra parameter k, 

h = N k (e- v2 / kv o - e -*~/**y (A3) 
with normalization factor 

N^ 1 = 4tt^ dvv 2 (V« 2 /^ _ e -t£./**o) (A4) 

which has no closed-form expression except for integral 
values of k. 



Appendix B: XENON100 acceptance 

In XENON100, the relative scintillation efficiency £ e ff 
is important for determining the mean number of SI pho- 
toelectrons corresponding to a given nuclear recoil en- 
ergy, through the relation 



(Bl) 



where S ee = 0.58 and S nr = 0.95 [79] are the suppression 
in the scintillation yield for electronic and nuclear recoil 
respectively, and the light yield L y = 2.20 PE/keV C0 is 
normalized at 122 keV ee [80] ■ The uncertainties in these 
parameters are negligible compared to that of £ c ff ■ 

In order to use the experiment to constrain DM mod- 
els, one could follow the likelihood method introduced in 
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|62j . A simpler method (similar to the one we adopt) is 
described in the thesis [61] , which requires the acceptance 
function A(E R ), that enters in the number of expected 
DM events in a given energy range (E\,E2) as 



N = £ 



E-2 



dR 
dEp 



A(E R )dE R 



(B2) 



where £ = 4843.2 kg-d is the exposure. The acceptance 
function takes account of £ c ff, the Poisson fluctuations 
of the 51 signal, and the acceptance of the SI and the 
nuclear recoil cut Ccut(rii) [23:, 



A(E R ) 



si max 

E 

n-i=Sln 



P(n 1)M si(^))CcutK) (B3) 



where fisi(E R ) is given in eq. (Bl| and the SI window 
used by XENON100 is Sl min = 4, Sl max = 30. The 
acceptance of the applied cuts ( cut is the product of the 
data quality cut and S2/S1 electron recoil discrimination 
cut, shown in fig. 2 of ref. [23 . The data quality cut 
depends upon the DM mass, and we interpolate between 
the curves given in [23] • P is the Poisson distribution, 



P (ni , msi (E R )) = [J>si (E R ) ni e~' lsl ^ /m ! 



(B4) 



To compute an upper limit from the number of events, 
we use the maximum gap method of ref. |63j . which as 
described in sect. |3.3| restricts N < 3.1 at 90% c.l. in the 
energy window 8.4-12.1 keVnr. 
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